Doping - dependent anisotropic superconducting gap in Nai_ ( 5(Fei_ a; Co x .) As 
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The London penetration depth was measured in single crystals of self-doped Nai-^FeAs (from 
under - to optimal - doping, T c from 14 to 27 K) and electron-doped Na(Fei_ a; Co a; )As with x ranging 
from undoped, x — to overdoped x = 0.1. In all samples, the low - temperature variation of the 
penetration depth exhibits a power-law dependence, AA(T) = AT 1 , with the exponent that varies 
in a dome - like fashion from n ~ 1.1 in the underdoped, reaching a maximum of n ~ 1.9 in the 
optimally doped and decreasing again ton~ 1.3 on the overdoped side. While the anisotropy of the 
gap structure follows universal dome-like evolution, the exponent at the optimal doping, n ~ 1.9, is 
lower than in other charge - doped Fe-based superconductors (FeSCs). The full - temperature range 
superfhud density, p 3 (T) = (A(0)/A(T)) 2 , at the optimal doping is also distinctly different from 
other charge - doped FeSCs but is similar to isovalently - substituted BaFe2(Asi_ a: P :E )2, believed to 
be a nodal pnictide at the optimal doping. These results suggest that the superconducting gap in 
Na(Fei_ :E Co :E )As is highly anisotropic even at the optimal doping. 



PACS numbers: 74.70.Xa,74.20.Rp,74.62.Dh 

The Fe-based superconductors (FeSCs) represent a rich 
variety of materials [IH5] . Despite this chemical diversity, 
a generic trend is that a superconducting "dome" in the 
temperature versus doping phase diagram forms in the 
vicinity or coexisting with the long range magnetic order. 
Instability of electronic system caused by the proximity of 
magnetism and superconductivity is considered as a key 
to understand pairing mechanism in the FeSCs 1 . Con- 
trary to the cuprate high-temperature superconductors, 
which show robust nodal d-wave superconducting gap for 
all doping types and levels, the gap structure of FeSCs 
is strongly doping-dependent varying from isotropic to 
nodal gap structure [3HH] . An important question at this 
moment is whether there is any universality in doping 
dependence of the gaps among different FeSC materials. 

The universal trend was suggested in studies of 
122 family. In particular, the gap structure of 
both electron-doped Ba(Fe 1 _ x T x ) 2 As 2 ("BaT122", 
T = transition metal) [6TfH2] and hole-doped 
(Bai_ x K a: )Fe2As2 ("BaK122") [HHH] are nodeless 
and isotropic at the optimal doping and become 
highly anisotropic (perhaps nodal) at the dome edges. 
The isovalent-doped LiFe(As,P) [T7J [TS] and layered 
Ca 10 (Pt 3 As 8 )[(Fei_ :c Pt :c ) 2 As2]5 ("CalO-3-8") com- 
pounds [JJ5] obey this trend as well. The only exception 
to this universal trend is found in isovalent - doped 
materials such as BaFe 2 (Asi_^P 2; )2 ("BaP122") [20"ll2T] 
and Ba(Fe 1 _ a ,Ru 2; )2As2 22J where nodal behavior was 
concluded even at the optimal doping. It was suggested 
by Kuroki et al. that nodal vs. nodeless behavior is 
controlled by the pnictogen height at least in a 1111 
family [23] and Hashimoto et al. extended this sugges- 
tion to 111 and 122 families [18]. It is also believed that 
the same parameter controls the transition temperature, 



T c [2~31 |2"5] . Other factors that might control nodal 
vs. nodeless gap are doping-dependent electronic band 
structure and competition between intraband and 
interband interactions, taking into account different 
orbital contents in each band (3] [4j [26]. Therefore, it is 
important to study other Fe-based systems to test these 
ideas and determine possible universal trends, if they 
exist. 

NaFeAs is one of few stoichiometric FeSCs with T c ~ 
12 K 27J, isostructural with LiFeAs (T c ~ 18 K). Ap- 
plication of 3 GPa pressure raises T c to 31 K, indicat- 
ing its effective-underdoped nature [3H] distinct from the 
slightly overdoped nature of stoichiometric LiFeAs |29j . 
Supporting this underdoped nature of parent NaFeAs, 
coexistence of magnetism and superconductivity was 
found in NMR, neutron scattering, electronic transport 
and ^(SR measurements [30453] . The coexistence region 
extends to optimal-doping in NaFeAs with electron Co- 
substitution Fc, 31 . The angle-resolved photoemission 
spectroscopy (ARPES) study in near optimally-doped 
NaFeo.95Coo.05As revealed the quasi-nested Fermi sur- 
faces connected by Q = (jr, it) wave vector, suggesting an 
important role of magnetic instability in superconductiv- 
ity pH], similar to the 122 family. This is in stark contrast 
to the absence of nesting in LiFeAs [35], despite recent 
reports from inelastic neutron scattering, indicating in- 
commensurate magnetic correlation in LiFeAs [36 . The 
gaps obtained in the ARPES study in NaFeo.g5Coo.05As 
are isotropic [33]. However, they persist well above T c 
and thus are unlikely related to superconductivity. Scan- 
ning tunnelling spectroscopy revealed V - shaped curves 
[37j , which may be consistent with both nodal and node- 
less gaps. Specific heat measurements revealed a jump 
at T c , but did not extend to low enough temperature 



2 



to reveal characteristic nodal vs. nodeless behavior [55]. 
Considering the ambiguity of the phonon contribution 
subtraction and multi-gap nature of superconductivity in 
this compound, in the current stage the specific heat re- 
sults cannot rule out the existence of (at least one) highly 
anisotropic gap. Therefore the important question on the 
anisotropy of superconducting gaps in NaFei-^Co^As re- 
mains open. 

In this work, we study doping-cvolution of the super- 
conducting gap structure in NaFeAs single crystals us- 
ing a high resolution tunnel-diode resonator (TDR) tech- 
nique. The doping level was controlled with either self- 
doping in Nai-^FeAs caused by oxidative environmental 
deintercalation (underdoped to optimally doped range) 
or Co-doping in Na(Fei_ x Co a; )As (spanning the whole 
superconducting dome). We found that, similar to most 
FeSCs, the variation of the London penetration depth is 
best described by a power-law function, A A = AT n , with 
the exponent n and pre-factor A varying systematically 
with doping in a dome-like fashion. The values of n in the 
underdoped (n ~ 1.1) and heavily overdoped (n ~ 1.3) 
compositions strongly suggest the superconducting gap 
with line nodes. Even at the optimal doping the expo- 
nent n ~ 1.9 is notably lower than in the other charge 
- doped FeSCs [5L The superfluid density, p s (T), at the 
optimal doping is also incompatible with isotropic full 
gap. Instead, p s {T) at x = 0.05 is remarkably similar 
to that of optimally isovalent-doped BaFe2(Asi_ ;r P :r )2, 
a known nodal pnictide superconductor even at the opti- 
mal doping [20] ■ Our observations suggest that, despite 
the universal tendency for highest anisotropy at the dome 
edges in LiFe(As,P) [T7|[TS], BaCol22 0, BaK122 [T3J 
and a layered CalO-3-8 system [19], the gap structure 
at the dome center can highly anisotropic even in the 
charge-doped systems. 

Single crystals of Na(Fei_ a; Co a: )As with x = 0, 0.02, 
0.025, 0.05, 0.08, and 0.10 were synthesized by sealing in 
the mixture of Na, Fe, As and Co together in Ta tubes 
and heating at 950 °C, followed by 5 °C/hour cooling 
down to 900 °C [33] ■ We determined x as used in load 
to crystal growth. The samples were stored and trans- 
ported in sealed containers filled with an inert gas. Sam- 
ple preparation for TDR measurements was done quickly 
in air within about 5 minutes to minimize uncontrolled 
environmental exposure which can induce an increase of 
T c @D1[H]. The process of bulk "self-doping" can be well- 
controlled using reaction with Apiezon N-grease, and can 
be facilitated by ultrasonic treatment stimulating dein- 
tercalation of Na + ions. Duration of the sonication treat- 
ment was used to control T c of sample, see Ref. [JT] for 
details. The low-temperature variation of the in-plane 
London penetration depth, AA(T), was measured using 
the TDR technique described elsewhere 02] [43] . 

Figure [I] shows the low temperature variation of AA(T) 
in two samples A and B of Nai-^FeAs crystals, which 
shows a clear evolution from almost T-linear (pristine) to 
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FIG. 1. (Color online) Low temperature part of AX(T/T C ) 
in two Nai-^FeAs samples (A and B) in pristine and ultra- 
sonic treated (one-hour) states. Dashed lines are representa- 
tive power-law fits up to T/T c = 0.1 for quantitative analysis. 
Inset shows full-temperature range variation of AA(T) for the 
same samples. 




FIG. 2. (Color online) (a) Low-temperature part of AX(T/T C ) 
in NatFei-zCo^As for x = 0.02, 0.025, 0.05, 0.08, and 0.10. 
Dashed lines are representative power-law fits up to T/T c = 
0.2 for quantitative analysis. Fitting results are summarized 
in Fig. [3] Inset shows full-temperature range variation of 
AA(T) for the same samples. 

almost quadratic (one-hour ultrasonic treatment) behav- 
ior. Measurements were done immediately after opening 
the ampoule (pristine state) and repeated after one-hour 
ultrasonic treatment, which increased T c to 22 K (See 
the inset). Successive treatments increased T c to a max- 
imum value of ~ 27 K (after 2 and 3 hour treatments), 
but longer treatments lead to sample degradation. 

Figure [2] shows the low temperature variation of AA(T) 
in Na(Fei_ x Co 2; )As for x = 0.02, 0.025, 0.05, 0.08, and 
0.1. Judging by T c (x) in the inset, the samples cover 
the superconducting dome from lightly underdoped to 
heavily overdoped regime. For a quantitative analysis of 
low temperature behavior, a fit to power-law function, 
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FIG. 3. (Color online) Phase diagram of NaFeAs system, 
(a) superconducting "dome", T c (x), (b) power-law exponent 
n(x), and (c) prefactor A(x). Self-doped results were placed 
against sonication time (top-axis) and scaled by their T c (see 
text). Circles show the data for Co-doped samples. Other 
symbols show the data for self-doped samples. 



AA(T) = AT n , was performed over a temperature inter- 
val from the base temperature up to several upper limits, 
varying from T c /10 to T c /5. The representative fits for 
T c /5 limit are shown as dashed lines. 

Figure [3] presents a summary of our results for NaFeAs 
system. A doping dependence of T c (panel a) in 
Na(Fei_ x Co x )As is shown against actual x (bottom 
axis). The data for Nai-^FeAs are plotted against soni- 
cation time (top axis) after scaling that the maximum T c 
is matched with that from the Co-doped case (i.e., 2.5- 
hour sonication treatment equivalent to x — 0.025). The 
parameters of power-law fit, AA(T) = AT n , are summa- 
rized in Fig. [3] (b) and (c). For each data set, fitting was 
performed for several temperature ranges to examine the 
robustness of the fitting parameters. In Nai-^FeAs, the 
exponent n is about 1.1 in the most underdoped com- 
position (pristine samples) and increases to 1.8 after 3 
hour sonication when the highest T c 27 K is reached. 
Similarly, larger exponent, n ~ 1.9, was obtained in 
Na(Fei_ 2; Co a ;)As for x = 0.025 (optimal doping) and 




FIG. 4. (Color online) Experimental superfluid densities, 
p s {T), of optimally-doped FeSCs. Top four curves for LiFeAs 
[17] , BaCol22 [5], BaK122 [13], and CalO-3-8 [19] were off- 
set for clarity. Red solid curves show the best fits to an 
isotropic s-wave two-gap 7-model [J3] demonstrating good 
agreement with the data. In a clear contrast, the p s for op- 
timally - doped NaFeo.975Coo.o25As (green open-squares) and 
NaFeo.95Coo.05As (blue open-circles) cannot be fitted with 
this fully gapped model, especially at the lowest tempera- 
tures, zoomed at in the inset. For comparison, we also show 
the data for optimally - doped BaP122 (x = 0.33, black solid 
dots) [5D] as well as single d— and s— gaps (dashed lines) in 
both clean and dirty limits. 



it monotonically decreased to n ~ 1.3 for overdoped 
x = 0.10. In general, the results for both Nai-^FeAs 
and Na(Fei_a;Co x )As follow universal trend that the gap 
anisotropies are highest at the dome edges and smallest 
at the dome center. The values at the dome edges are no- 
tably below the lowest possible value (n ~ 1.7) expected 
from s± scenario with strong pair-breaking scattering [5] , 
suggesting the formation of nodes. In addition, the low 
value of exponent n ~ 1.9 at the dome center makes a 
clear distinction of Nai-^Fei-xCo^As from the other 
charge - doped FeSCs. The exponents close to n ~ 2 
can be explained by strong pair breaking effect in s± or 
ci-wave pairing scenarios. 

For a more complete analysis of superconducting gap 
structure, especially in a multi-gap system, the low- 
temperature power-law behavior is not sufficient, as 
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it only reflects quasiparticle excitations around gap 
minima. Figure [4] shows normalized superfluid den- 
sity, p s (T) = (A(0)/A(T)) 2 , plotted against T/T c 
for several optimally doped FeSCs. To calculate 
p s in NaFeo.975Coo.o25As (green open-squares) and 
NaFeo.95Coo.05As (blue open-circles), we used the val- 
ues of Ao = 410 nm and 354 nm, respectively, obtained 
from /xSR measurements |31j . The p s for other mate- 
rials (offset vertically for clarity in Figure [4]) include 
LiFeAs Q2], BaCol22 0, BaK122 [T3], and CalO-3-8 
[TjJ], all measured in our group. Red solid curves are 
the best fits to a self-consistent isotropic s-wave two- 
gap 7-model [44] showing good fit quality in all four 
cases. In stark contrast, an attempt to fit the data for 
NaFe .975Co .o25As and NaFeo.95Coo.05As (both are near 
the optimal doping) with the same 7— model clearly fails 
especially at the low-temperature end, where the distinc- 
tion between full and nodal gaps is most pronounced, as 
zoomed in the inset of Fig. |4j For comparison, we also 
show p s (T/T c ) expected for a single-gap s- and d-wave 
pairings both in clean and dirty limits (dashed curves). 
These single-gap curves clearly deviate in the whole tem- 
perature range, confirming multi-gap nature of supercon- 
ductivity in Na(Fei_ a; Co a ;)As (as well as other FeSCs). 
Surprisingly, p s {T/T c ) of NaFeo.95Coo.05As is virtually 
the same as that of optimally - doped BaP122 (x = 0.33, 
black solid dots), a known nodal isovalent substituted 
pnictide |20j . and shows notable difference from the ex- 
ponential saturation in a fully gapped LiFeAs [TT1 [18] . 

Our results show that the superconducting gap even 
at the optimal doping is highly anisotropic in charge- 
doped Na(Fei_ a; Co a; )As. Similar to other FeSCs it be- 
comes more anisotropic upon departure from the opti- 
mal doping. To explain such sensitivity to the doping 
level and appearance of nodes in some materials, several 
proposals have been put forward, (i) Following a sug- 
gestion by Kuroki et al. that nodal vs. nodeless behav- 
ior is controlled by pnictogen height in 1111 family [25] , 
Hashimoto et al. extended this suggestion to 111 and 
122 families [TS]. (ii) Nodal structure of the supercon- 
ducting gap was linked with the evolution of the Fermi 
surface topology and orbital content in different bands 
H [2S EE! ES]. The bandstructure of NaFeAs is signif- 
icantly different from 122 compounds, which itself varies 
significantly between different hole, electron and isovalent 
substitutions with no apparent correlation that may ex- 
plain appearance of nodes only in some compounds, (iii) 
Another approach is to incorporate changes in both, the 
bandstructure and renormalization of intra-band and in- 
band Coulomb interactions. According to Ref. 5] (1) for 
weakly and moderately electron-doping, the propensity 
of s± and d— wave pairing channels are comparable; (2) 
weakly and moderately hole-doping, s± pairing is dom- 
inant; and (3) for strongly electron- and hole-doping, 
d-wave channel is dominant. From this classification, 
regime 1 seems to be most relevant for the optimally- 



doped Na(Fei_ a; Co a ;)As and regime 3 is appropriate for 
the overdoped compositions, consistent with our results. 

In conclusion, a combined study of self-doped 
Nai-^FeAs and charge-doped Na(Fei_ 3; Co a ;)As estab- 
lishes the systematics of the superconducting gap evo- 
lution with doping in the 111 compounds. Similar to 
other FeSCs, we find a universal tendency of developing 
highest anisotropy at the dome edges. However, we show 
that, even at the optimal doping, the low temperature 
penetration depth study indicates highly anisotropic su- 
perconducting gap. In addition, temperature - dependent 
superfluid density at the optimal doping is remarkably 
similar to p s (T) of isovalent-doped BaFe2(Asi_ ;E P ;E )2, a 
known nodal pnictide superconductor at the optimal dop- 
ing. 

After posting the preprint of this paper [57], two 
studies using crystals from Prof. X. H. Chen group 
[37l 138] suggested full isotropic superconducting gap in 
NaFei-^Co^As at optimal Co-doping. To check if the 
difference in the conclusions is caused by the difference 
in the samples of this reactive material, we asked Prof. 
X. H. Chen to send us crystals for a cross - examination. 
We received and measured optimally doped samples with 
x=0.028. The results were virtually identical (within the 
noise level) to those presented in Fig. [3] and |2J ruling out 
the sample dependence of our conclusions. Possible rea- 
sons for the discrepancy with other studies are discussed 
in the introduction. 
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